function [ output ] = Gauss( x,ampl,center,sigma )
 output= ampl*exp(-((x - center).^2)/(2 * sigma^2)); 
end

